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Abstract 

The study of the two shell system started in our first paper “Pair of null 
gravitating shells I” is continued. An action functional for a single shell due to 
Louko, Whiting and Friedman is generalized to give appropriate equations of 
motion for two and, in fact, any number of spherically symmetric null shells, 
including the cases when the shells intersect. In order to find the symplec- 
tic structure for the space of solutions described in paper I, the pull back to 
the constraint surface of the Liouville form determined by the action is trans¬ 
formed into new variables. They consist of Dirac observables, embeddings and 
embedding momenta (the so-called Kuchaf decomposition). The calculation 
includes the integration of a set of coupled partial differential equations. A 
general method of solving the equations is worked out. 


1 Introduction 


The present paper is the second in a series dedicated to the two-shell system. The 
first paper, Ref. [Q, will be referred to as I henceforth. The general motivation 
and aims are explained in I. Some solutions to the classical equations of motion 
containing two shells have been studied in I with two main results: First, certain 
coordinates have been chosen on the space of all solutions; they are candidates for 
a complete set of Dirac observables of the system. Second, all symmetries of the 
system that are associated with diffeomorphisms between solution spacetimes have 
been found. These symmetries will dehne the most interesting observables as well 
as the true Hamiltonian of the system. 

To proceed, the space of solutions ought to be promoted to a physical phase space, 
which is a symplectic manifold. Without knowledge of the symplectic structure, not 
even the decision can be made whether some physical coordinates are differentiable. 
Therefore, we turn in this paper to the canonical theory. 

As far as we know, no action functional for the two-shell system exists in lit¬ 
erature. For a single spherically symmetric null shell, a Hamilton action principle 
has been studied by Kraus and Wilczek 0 and by Louko, Whiting and Friedman 
ip. We shall start from this action and generalize it for two and more shells. The 
generalization is easy, and a proof is given that the resulting equations of motion 
admit the correct solutions for the system, including all solutions^ described in I. 
This, of course, is not a full justification for the choice of action. As it is well known, 
given dynamical equations admit many non-equivalent action principles from which 
they follow. An example for a single massive shell is investigated in [Q. A direct 
reduction of the Einstein-Hilbert plus a shell-matter action to spherical symmetry 
might lead to a non-equivalent theory. Ref. is an example of such an action prin¬ 
ciple without symmetry, but it is valid only for a single shell. The generalization of 
the Louko-Whiting-Friedman action to more shells and intersections is simpler than 
that of the action in so we have chosen to start from it. We assume that both 
actions lead to equivalent theories. 

The plan of the paper is as follows. In Sec. we propose our generalization of the 
Louko-Whiting-Friedman action principle to include two and, in fact, any number 
of in-going and out-going shells. We prove that the variation of this action gives the 
equations of motion that we need. 

Sec. 1^ begins with the calculation of Poisson brackets between the parameters 
that have been chosen as coordinates on the space of solutions in 1. The strategy is 
to pull back the Liouville form defined by our action to the constraint surface, and 

^ These solutions are a subset of the solutions admitted by the equations of motion due to the 
additional assumption of the regular center on the left. 
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to transform variables so that the new ones consist of the parameters, embeddings 
and embedding momenta. This is the so-called Knchaf decomposition stndied in P . 
A similar transformation has been accomplished in []^ for the case of a single shell. 
The calcnlation inclndes integration of a set of partial differential eqnations of the 
hrst order. The integration method nsed in [0 seems to be rather closely associated 
with the particular problem studied and the gauge chosen there. We shall hnd a 
general method that works for any number of shells and in any double-null gauge. 


2 Canonical formalism 


2.1 The action 


The action functional for a single light-like shell and its gravitational held has been 
written down in Ref. |^. It is modihed below so that it describes two such shells: 


& = dr 


pif 1 -h P2r2 + / dp{PAA + PrR - H 2 


( 1 ) 


The Hamiltonian has the same overall form as in 


H2=AfH + J\fpnp + Afo,E^ , 


but the new constraints are 


n 


2R2 

hiPi 

A 

PrR'- 


PaPr ^ RR" 


R 


A 


RR'K' ^ R'2 
A2 ^ ^ 


5(p- ri) + 


p;a 


A 

pi(5(p- ri) 


1 - 2 ) , 
P2^(p 


A 

1 - 2 ) • 


( 2 ) 

( 3 ) 


The momenta of the shells are pi and p 2 , their radial coordinates are ri and r 2 and 
their radial directions are rp and 772 - The dot denotes the derivative with respect to 
r and the prime that with respect to p. The “volume” variables A, R, Pa, Pr, M 
and are the same as in [Q and 0]. The meaning of the variables A, R, M and 
can be inferred from the spacetime metric 

ds^ = -M^dr^ + k^{dp + U^drY + R^dn^. ( 4 ) 


The momenta conjugate to the conhguration variables A and R are 


Pa =- WR') , 

and 

Ph = - A/-PP') - |.|A - (A/'-’A)'! . 


( 5 ) 

( 6 ) 
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Compared to ref. p], the only change is the presence of one more term in the 
Liouville form and of a corresponding term in each constraint. Each of these terms 
has the same form as in Q, but now, they carry indices so that the two shells can 
be distinguished. 

It is important to state the differentiability conditions on the shells, especially at 
a possible shell crossing (if r]i ^ 772 )- One can assume as in Q that the gravitational 
variables are smooth functions of p, with the exception that A/”', (A/"^)', A', i?'. Pa 
and Pr may have hnite discontinuities at isolated values of p. The coordinate loci 
of the discontinuities are smooth functions of r for each shell. This follows from (i) 
the conditions at isolated single shell points which are the same as in |^, (ii) from 
the conditions at any shell crossing as described in Sec. 2.2 of I and (iii) from the 
corresponding choice of foliation: the metric with respect to coordinates r and p 
may be piecewise smooth and everywhere continuous. 

2.2 Equations of motion 

The rest of the section will be devoted to a check that the variation of the action ( 0 ) 
yields the proper dynamical equations for two shells. These equations have to imply 
(a) that the geometry outside the shells is the Schwarzschild one and (b) that the 
shells move along null surfaces. These two properties, together with the continuity 
of the metric everywhere and the assumption of a regular center on the left, have 
been used in Sec. 2 of I to construct all solutions. The continuity of the metric is 
part of the dehnition of the conhguration space, and the existence of a left regular 
center is an assumption of our model, but the two other properties have to result 
from the variation of the action. 

The continuity of the metric dehnes a class of coordinates; we shall work with 
such coordinates in this section. We also assume that the foliation coordinates r 
and p are C^. The variations of the action (|l]) with respect to the gravitational 
variables A/", A/"^, P, A, Pr and Pa give the constraints, 

= 0 , Pp = 0 (7) 
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and the dynamical equations 


( 8 ) 

(9) 


A 

R 

Pa 


Pr 


M 


APa Pr 


R 

ATPa 




R 

2 

M'RR' 


Aff’R' , 

-W~[a) + + 


2?72P2 

A2 


5(p-r2) 


+ Af^Pk , 


Af 


A2 

API P^P, 


R^ 


A-Tfi 

R^ 




'APR" 


+ {U^Pr)' . 


( 10 ) 

( 11 ) 


Outside the shells (p 7 ^ ri,p 7 ^ r 2 ), they coincide with the equations of [Q and 
imply that the geometry there is the Schwarzschild one corresponding to some, as 
yet arbitrary, value of Schwarzschild mass parameter. 


2.2.1 Outside the crossing 

Let us hrst consider a Cauchy surface that does not contain a crossing point of the 
two shells, ri 7 ^ r 2 . This includes data for crossing as well as for parallel shells. 

The equations that we obtain at p = r*, i = 1,2, result partly from setting the 
coefficients at the delta functions (5(p — ri) and 5(p — r 2 ) in Eqs. 0-(Q) equal to 
zero, and partly from the variations of the action (|lD with respect to rj and p*. If 
we use the notation 5r(-^) for the coefficient at the 5{p — r) in the expression X, 
then the equations read: 


0 

0 

5r.(PA) 

^rAPR) 

i-i 

P^ 



+ 


ViPi 

Ap=ri 




5ri{APk) + Pi , 



-iAf)p=r. 


ViPi + (AA)p=rAriiPA) ) 




+6rA[Af^PRr) , 


( 12 ) 

(13) 

(14) 

(15) 

(16) 
(17) 


By inspection, the equations with any hxed value of i coincide with the one-shell 
equations of |^. Thus, outside the shells, as well as at the shells outside the crossing. 
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we obtain equations that have been shown in to imply the required properties 
(a) and (b). 

To obtain expressions for the momenta pj conjugate to r^, we can use the con¬ 
straints (1^) and (|T^) because we shall always work only at the constraint surface 
within the present series of papers. We obtain the equations 


R{ri)Ari{R')+r]iPi = 0, (18) 

A(ri)Ari(PA) + Pi = 0. (19) 

Here, we use the notation Aa.(X) for the jump limp^ 3 ._|_ X — limp^ 3 ._ X at the shell 
point p = X. 

In any double-null coordinates U and V, the metric has the form 


ds^ = -AdUdV + R^d^l^ 
and the transformation formulae are: 

A = V-AU'V , 

iiv - vu' 


M = 
MP = 


2U'V 

UV + VU' 
2U'V' 




It follows that 


( 20 ) 


( 21 ) 

( 22 ) 

(23) 


E! = RjjU' + RyV' , 


Pa = - RvV) . 

The last equation is derived from Eq. (|^) and 
following simple Lemma. 


(21)-(23 


(24) 

(25) 

We shall often need the 


Lemma 1 Let (p be a continuous function of U and V and let its derivatives have 
jumps only at U = U{r) and V = H(r), being themselves continuous elsewhere. 
Then ip^u is a continuous function of V everywhere except for U = U{r) where the 
jump A{ip^u) is a continuous function ofV for all V (even atV = V{y)). Similarly 
for ify. 

The Lemma follows immediately from the continuity of p>. 

Applying the Lemma to the function R, which is continuous, we find that A{R ,u) 7^ 
0, A{Ry) = 0 for Pi = -|-1 and A{Ry) = 0, A(Ry) ^ 0 for pi = —1. 

Substituting this into Eqs. ([T8|) and (|^, one finds the desired expressions: 


Pi = -R{vi)A^.{Ry)U' 


(26) 
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for rji = +1 and 


p, = i?(r,)Ar,(i?,v)f/' 


(27) 


for rji = —1. The orientation of the functions U, V, r and p is such that U' and 
V > Q (see [0]). 

2.2.2 At the crossing: the constraints 

Let us now turn to the Cauchy surfaces that contain a crossing point of the shells, 
where we have ri = r 2 = r. Here, there is only one ^-function, (5(p — r), and there 
is only one coefficient in all variation equations at — r). Thus, instead of two 
equations, we obtain only one in each case. A new analysis is in order. 

Let us first consider the constraints (i- The coefficients of S{p — r) in both 
constraints are 

m 


W) = m + + 


K{R') , 


A(r) A(r) A(r) 

6r{n^) = -Pl - P 2 - A(r)Ar(PA) • 

To study the equations SriH) = 0 and SriH^) = 0, it is advantageous to choose 
some double-null coordinates in which the metric is continuous. Similar tactic has 
been pursued in the Appendix B of 0. The expressions (|^) and (p5[) for R' and 
Pa in terms of general double-null coordinates yield now 

Ar(P') = Ar(Pc/)P'(r) + Ar(Py)H'(r) , 

P(r) 


Ar(PA) = 


A(r; 


[Ar(Pf/)P'(r)-Ar(Py)H'(r)] . 


Thus, the equation SriTi.) = 0 becomes 

P(r)Ar(P,[/)P'(r) -h P(r)Ar(P,y)H'(r) = -pipi - P2P2 
and = 0 is 

P(r)Ar(Pf/)P'(r) - P(r)Ar(Py)H'(r) = -pi - p2 . 

Let Pl = —p 2 = 1- Then, we obtain from Eqs. (^8]) and (|^) 

Pl = -P(r)Ar(Pc/)P'(r) , 

P2 = P(r)Ar(Py)H'(r) . 

If Pl = —p 2 = — 1, then 

= P(r)Ar(Py)H'(r) , 

= -P(r)Ar(Pc7)P'(r) . 


(28) 


(29) 


Pl 

P2 
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This can be summarized by 


Pout = -R{r)Ay,{Ru)U'{r) , 
Pin = Rir)AriRy)V'ir) . 


(30) 

(31) 


Notice that these equations are exactly the same as Eqs. (^) and (0) for shells 
outside the crossing. They can serve as dehnition of the momenta everywhere. From 
Lemma 1 and Eqs. (^) and (|3l|) we can infer that Pout is continuous along U = U{v) 
and Pin along V = V(r). 

The fact that we obtain only two relations from the four equations 5^(7^) = 0, 
Sr 2 iR) = 0, Sr^iTLp) = 0 and — 0 outside the crossing is due to the mutual 

dependence of Ti- and Tip-equations (see 0) in this case. At the crossing the Ti- 
and Tip-equations are independent. 


2.2.3 At the crossing: the shell variables 

Next, we study the variation of the action (|^ with respect to the shell variables r* 
and Pi. The variation is to be done so that ri and r 2 are considered as independent 
variables even where ri = r 2 ; therefore, one must vary Erst and only then set ri = r 2 . 
The resulting Euler-Lagrange equations then are 


M{t) 

= TTT^* 

A(r) 




and 


Pi = 


- ( 


\Hp) 


Pi 


(32) 


(33) 


p=r 


Eq. (|3^ is analogous to Eq. ([I6|) and Eq. (|33|) to Eq. ([TtI) . 

Let us assume that the index i is chosen in such a way that rji remains constant 
along the trajectory of the Ath shell for each i = 1, 2. This is not the only possible 
assumption if the shells cross each other. For example, one could choose i = 1 for 
the innermost shell and i = 2 for the outermost one. Then the subsequent argument 
has to be modified. 

The interpretation of Eq. (|3^) is the same as in |^, if we use our convention on the 
index i. The right-hand side of Eq. (^) is then continuous along each shell because 
the functions A/”, A/”^, and A are components of the metric in the coordinates r and 
p. The vector tangential to the shell, 

/f = (l,i-„0,0), 

is a null vector as a consequence of Eq. (^). Hence, the equation completes the 
dynamical information that is necessary for construction of solutions. The remaining 
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equations for the shell have, therefore, to follow from Einstein’s equations outside 
the shells, Eqs. (|30|), (|^ ) and 

Let us extend the dehnition of the quantity /f to a whole neighborhood of i-th 
shell by 

V(P) ■ 


The function /f(p) can be calculated in terms of the metric and embeddings for an 
arbitrary double-null gauge. Using the transformation Eqs. (0)-(ED. we obtain 


;f(p) = -L (34) 

for rji = +1 and 

‘f(p) = ~ (35) 


for rji = —1. At the shells, this follows also directly from the relations t/(r, r(r)) = 
const for rji = -|-1 and V (r, r(r)) = const for rji = —1, and it is analogous to Eq. (B3) 
in 0. The right-hand sides of Eqs. (^41) and (pSD are, of course, continuous functions 
in a neighborhood of Uth shell because we have required the foliation to be C^. 

Let us turn to Eq. (|33|) . As in [^, the hrst question is whether the right-hand 
side is continuous across the shells; if it is not we have a bad ambiguity. Let rji = +1 
and let us consider the question of derivatives of the function /f(p) with respect to 
p from both sides of Ath shell (the argument is completely analogous to that in 
Appendix B of [|]): 


We have, however. 



1 

U' 



p=r± 



p=rit 


(f>y.r± + 

Av(r,r(r)). 


and this is, of course, well dehned. Hence, Eq. 

= — 1 is similar. 

Next, we show that Eq. (|3^ is satished if Eqs 
R is continuous at the crossing 


is not ambiguous. The proof for 


J) and (1^ hold and the function 
As we have seen. Lemma 1 then implies that 
the jumps A(/2 f/) and A{Ry) are continuous along the corresponding shells, even 





through the crossing point. Let = +1 and let us substitute Eq. (p^ ) for p* into 
Eq. (|33[). A simple calculation gives 

L|A,((fi=),t,)r(r,r(r))=0 

along the shell and also at the crossing point. However, outside the crossing, this 
equation follows from other equations (see [Q), so it does not contain any new 
information, even at the crossing: if Ar((i?^),f/) is time independent from both sides 
and continuous everywhere, then it must be also constant everywhere. 


2.2.4 At the crossing: A and R 

The remaining dynamical equations for the shells are obtained by varying the action 
(P with respect to the variables A and R and setting the coefficients at — r) in 
the resulting equations equal to zero. The variations with respect to Pa and Pr give 
just dehnitions of Pa and Pr, which do not contain 5{p — r). Thus, at the crossing 
point, we obtain the equations: 

= ^^(hiPi + h2P2) , (36) 


and 

sAPr) - m-rMPk) + = 0 ■ ( 37 ) 

while there are four equations outside the crossing, (0 and (0), that can be cast 
as follows: 

. (38) 

Apr*) 


and 


■Sr,)/’*) - + ^^r.(V") = 0 


A(ri) 


A(ri) 


(39) 


for each i = 1,2. 

Let us study the four discontinuous functions R', Af', Pa and Pr, the r and p 
derivatives of which feature under the 6r signs in Eqs. (p6D-(p9|). It is advantageous 
to express them as functions of U and V and of the foliation U{t,p) and V{t,p). 
For P' and Pa, Eqs. (|2^) and (p^ yield the desired expressions. For Af, we have 
immediately 


Af = Af,uU' + AfyV . 


(40) 


9 








For Pr, Eq. (|]) and (pT[)-(p3[) yield 


- R,V' 


We observe that 

U" V" 
TF~V' 


In 




V' 


2A 


In 


2A 


2U' 2V' 




v 


u' + 


,u 


In 




V' 


V', 


,v 


where \a{—U'/V) is considered as a (continuous) function of U and V. Thus, we 
obtain, finally: 


Pr = 


2AR 


(r^a) u' —5— (r^a) V' 

V Ru 2AR^ Rv 


R 


In 


'-V' 

v 


J ,u 




In 


-f/'' 

V 


V' . (41) 


,v 


Eqs. (1^), (^51) , (]40|) and (41) show that each discontinuous function is a sum of 
terms of the form or ipRy, where (p and R are continuous functions of U and 
V. The jump structure of such terms is given by the Lemma 1. 

Let us denote by -Rout) -^out; -RAout and Pnout the sum of all terms in the right- 
hand sides of Eqs. (^), (|25|), (^0]) and (41) that contain only the t/-derivatives, and 
similarly for R{^ etc. Thus, e.g.. 


-RRout 


2AR 



[ y-u'\i 

(r^a) u' + - 

In T., 

V Ru 2 

V ; 


UR 


,u 


etc. Eqs. (pq) and (^) can then be written as follows: 


[Sr{{PAouty)-AfRr)Sr{{PAonty)] + [5r((PAin)-) " Ar"(r)5r((PAin)')] 

A/’(r 


A2(r) 


(houtPout + VinPin) , (42) 


and 


SAiPR^uA) - MymPRout)')) + RRsAiRLtY) + 

+ SAiPRiny) - mr)SA{PmnY) + ^A((KY) + = 0 , ( 43 ) 


while Eqs. (|38[) and (^) become, for the out-going shell, 

l 5 r...((-PAo„.)') - V'-(r„„.) 4 „,((FA„..)') = ^ 4 ^>(o„.Po„. , 


( 44 ) 
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and 


^rout ( (-P-Rout) ) (I’out) rout ((Pr out } 

M (^out) 

A (rout) 




^((A/'™.)') = 0 (45) 


and similarly for the in-going shell. The reason is that only the derivatives of the 
out-terms give contributions to 5’s along the out-going shell because the in-terms 
are continuous, according to Lemma 1. 

The left-hand sides of Eqs. and (45) are obtained from the jumps of IJ 
derivatives of continuous functions in a way that is continuous along the out-going 
shell. The jumps themselves are continuous along the shell because of Lemma 1. 
However, the out-parts of Eqs. (42) and (43) are made in the same way from the 
jumps of U derivatives at rout = r. Because of the continuity of all terms along 
the shell, the out-part of the left-hand side of Eq. (42) is the rout —^ r limit of the 
left-hand side of Eq. (|4^ . Analogous claims hold for left-hand sides of Eqs. (43) 
and (45) as well as for the in-terms. Moreover, because of the continuity of pout and 
r/out along the out-going shell, the out-term on the right-hand side of Eq. (42) is the 
limit of the right-hand side of Eq. (^4]). 

It follows that each of Eqs. (42) and (43) can be considered as sum of two equa¬ 
tions, one being the limit rout —r of the corresponding out-going shell Eqs. (@) 
and (45), the other being the same limit of an in-going shell equation. Hence, if the 
out- and in-going shell equations hold for all values of rout and rin outside the cross¬ 
ing point, then Eqs. (42) and (43) are also valid. This is implied by the continuity 
conditions on the phase-space variables. 

However, Eqs. (^41) and (45) are satished because they follow from Eqs. (^) and 
(p9|), and these, in turn, follow from other dynamical equations. This has been shown 
in for a single shell and the proof is, formally, the same in our case because all 
one-shell equations outside the shell crossing are formally identical with equations 

in g. 

To summarize: We have shown that the action (P gives proper dynamical equa¬ 
tions for the two-shell system. 


3 The Liouville form at the constraint surface 

Our hnal aim is to calculate the Poisson brackets between Dirac observables such 
as Mm, Mr, Vm 2 — and Vr defined in 1. Our method will employ the property of 
the pull-back Op of the Liouville form 0 to the constraint surface T that it depends 
only on the Dirac observables. Its external differential then defines the symplectic 
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form of the physical phase space. In the present section, we develop some general 
tools in this line. 

The Liouville form of the action (|^ can be written as follows: 

roo 

0 = Pifi + P2r2 - MooEoo + dp (PaA + PrR) ■ (46) 

Jo 

We have inclnded the bonndary part of the Hamiltonian into 0; the form of this 
part justifies such inclusion (cf. [§], [0). We can now start to transform (^) into the 
Kuchaf variables corresponding to arbitrary double-null gauge. While 0r does not 
depend of the gauge and dependent degrees of freedom, the old variables contain 
them and we must use the complete transformation; it goes from the variables R, A, 
Pr, Pa to observables, gauge variables and dependent variables. Still, the resulting 
0r contains then only the former and none of the latter. 

To write down such a transformation, we shall use a particular Kuchaf decom¬ 
position; for dehnition and existence, see @]. To this aim, we shall choose an arbi¬ 
trary double-null gauge, represented by the coordinates U and V. We also let the 
choice of Dirac observables open, denoting them by o^, fc = 1, 2,. .., 2N. The com¬ 
plete set of final variables consists, therefore, of the physical variables o^, the gauge 
variables represented by the embeddings {U{p),V{p)) and the dependent variables 
represented by the embedding momenta Pu{p) and Py(p); at the constraint surface, 
Pu{p) = Pv{p) = 0- 

To start the calculation, we just need to know that the metric (^) depends on 
oh 

A = A{U, K; o) , R = R{U, V; o) . 

Then we use the transformation formulae and the dehnitions and 

of the momenta. Such calculation has been already carried out in |0 so we just take 
over the relevant general formulae from there. 

The calculation in then leaves this general stage and proceeds by making a 
particular choice of o’s as well as of the gauge P and V. It grows rather complicated 
and can be accomplished by some miraculous tricks whose nature seems to be closely 
connected with the particular choices of gauge and observables made. The main 
purpose of the present section is to reveal a general structure that underlies the 
tricks and that is entirely general. 


3.1 The volume part 


The form (^) can be divided into a boundary part (the hrst three terms on the right- 
hand side) and the volume parts. Each volume part is associated with a particular 
component of the space between the shells; it has the form 

rb 


dp (TaA -|- PrR) , 


J a 
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where a and b are values of the coordinate p at the boundary of the volume. For 
example, a = 0 and 6 = ri, or a = r 2 , 6 = oo, etc. 

The only “volume variables” among the hnal set are the embedding ones, U{p) 
and V{p), and the embedding momenta Pu{p) and Pv{p)- However, ©r cannot 
contain U{p) and V{p) because they are gauge variables; still less can it contain 
Pu{p) and Pv{p) because they vanish at F. Hence, we expect that the volume parts 
of ©r can all be reduced to some boundary terms. A general account of such a 
reduction is now given. 

We can see that the form of the volume parts is independent of the system, more 
precisely, of the number of the shells in the system. Hence, we can use the methods 
of Ref. 0] for its transformation to the new variables. 

For any double-null gauge, the equations 


ARRyy -\- AR jjR y -f A 

= 0, 

(47) 

^P,uu ~ ^,uP,U 

= 0, 

(48) 

ARyy — AyRy 

= 0 

(49) 


represent the condition that the transformation is performed at the constraint sur¬ 
face F (cf. 0, Eqs. (32)-(34)). 

For the transformation of any volume part into the new variables, we can make 
the ansatz 

0alr= ["dpiifU + gV + h.dy + p] , (50) 

as in 0. In the expressions for the functions /, g, hi and p, we can separate the 
terms with p-derivatives of the embeddings U{p),V{p) from terms in which U{p) and 
V{p) are not differentiated. The form of the transformation between the variables 
A, R, P[y and Pr on one side and U, V, Pr, Py and o* on the other then imply that 

/ = +F(C/V.o‘) . (51) 

9 = PR\n(-^^+G{U,V,o') , ( 62 ) 

K = ^\n(-'h'^+H,(U,V,o‘) , (53) 

p = RR^uU'- RRyV-^{R^uU'+ RyV')\n 

-FU' - GV + (f){U, V, o^) , (54) 
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and (see Eqs. (45)-(48) of 0) 


Fy ~ = 

^^{2AR^uv — A,uRy — ^yR,u) , 

(66) 

Hiy — F^i = 

~-;^G^R,iU ~ ^,iR,U ~ ^,uR,i) ) 

(56) 

Hiy — G^i = 

— (2y4Fiy — A^iRy — AyR^i) , 

(67) 

0 = 

0 . 

(58) 


Eqs. (p5|)-(|58|) are valid for any double-null gauge. 

The ansatz (|50|) leads to the following transformation of the volume part: 

©air = ifU + gV + hid" - ^pb)p=k - {fU + gV + hid" - (pa)p=„ + ^ 

By ignoring the total time derivative we are thus left with an equivalent form that 
contains only boundary terms: 

©air = {fU + gV + Kd" - ipb)p=h - if if + gV + h,d" - . (59) 

The boundary values a and 6 of p depend, in general, on the time parameter r as 
they can also describe positions of shells. If the intermediate boundaries are chosen 
to coincide with the positions of shells, the total Liouville form 0 reduces to the sum 
of contributions from (i) p = (corresponding to points a, b at which the embedding 
U{p), V{p) intersects a shell), (ii) p = 0 (corresponding to the center i? = 0) and (hi) 
p —cx) (corresponding to asymptotic inhnity). The boundary between two adjacent 
spacetime regions A4k and M.k+i is either a light-like hypersurface (dehned by an 
in-going or an out-going shell) or a crossing point (dehned by the intersection of two 
shells of different p). 

In this way, ©r can be transformed to a sum of boundary terms, provided that a 
solution to the system of differential equations (^ ) - (|T^) can be found. Let us study 
these equations. 

3.2 Properties of functions F, G and Hi 

Let us hrst establish some general properties of the functions F, G and Hi that 
follow from Eqs. (|55D -(|^. 

3.2.1 Freedom in the functions F, G, Hi 

Eqs. (^51)-(^) form an inhomogeneous system of linear partial differential equations 
of hrst order. Eqs. (|47|) -(^9D constitute the integrability condition for the system. 
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The general solution of the system can be written as a sum of a particular solution 
and a general solution of the homogeneous system. Let F^, and be a particular 
solution, and let F^, and be a solution of the homogeneous equations: 



Fy-G]u =0, 

(60) 


o 

1 

(61) 


o 

1 

(62) 

Let us define the function Cij by 



Q, := HF - Hi 

(63) 

and study its properties. 

Immediately from the dehnition, we have 



Cij ~Cji , 

(64) 

and 

Cij,k T Cjk^i T Cfcjj 0 . 

(65) 


The derivative of Eq. ([6l|) with respect to o' with subsequent anti-symmetrization 
in the indices i and j yield 

Cij,u = 0 . 

Using Eq. (|62[) in a similar way gives 

Gj^v H ■ 

Hence, Cij depends only on Dirac’s observables. Eqs. (Q) and (^) imply then that 
there is a function Gi{o) such that 

Cij = Gij — Gj^i . 

Let us choose an arbitrary Cj(o). Then the functions F^, G^ and have to 
satisfy the equations 


Flv - G]u 

= 0, 

(66) 

Hlu-Fl 

= 0, 

(67) 

Hly-Gl 

= 0, 

(68) 

Tji . _ zyi , 

= Gij — Gj^i , 

(69) 


Again, this is an inhomogeneous linear differential system. The following is clearly 
a particular solution: 

= G^ = 0 , Hi =Gi . 
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Any solution of the corresponding homogeneous system has, however, the form 

F^ = W,u, G^ = Wy, H] = W,, 

where W is an arbitrary function of U, V, and oh 

We have shown: Let and be a solution to the system (|55| ) - (|^ . Then 

any other solution F, G, has the form 

F = F^ + W,u, G = G^ + Wy, Hi = H^ + W,i + Gi , (70) 

where W is an arbitrary function of the variables U, V and o*, and Q is an arbitrary 
function of Dirac’s observables. 

If we substitute the solution (|70|) into Eq. for 0a|r, we obtain for the terms 
containing the functions W and Cj (the sum of these terms is denoted by 50^ |r): 

50^|r = r dp [{WyU + WyV + WiO* + Qd^ + (-Wyl/' - W)-] . 

J a 

However, Dirac’s observables do not depend on p, so (Qd*)' = 0 and WyU' — 
WyV = tc', where we have dehned 

w { p , t ) := ID(f/(p,r),I/(p,r),o*(r)) . 


Using this, we have 

^©alr = / dp[{w)' - («;')■] = 0 

J a 

because all functions are G°° in the space between the shells. Hence, each solution 
( 0 ) leads to the same Liouville form. 

The W- and Cj-part of the hrst parenthesis in Eq. (|5^) is, however, a total time 
derivative: 


{WyU + WyV + + GV + Wyl/'i) - WyVb)p=b = 

dw{h,T)- dw{h,T) 


db 


-b + 


dr 


= ^«-(Kr).r). 


which, in general, is non zero. Similar result holds for the second parenthesis. Hence, 
different solutions lead to equivalent boundary Liouville forms. 


3.2.2 Gauge transformation of F, G and Hi 

Two different double-null gauges lead to two different transformation of a given 
volume part of the Liouville form. We can, therefore, ask how the functions F, G 
and Hi are transformed if the gauge changes. The leading idea of course is that the 
Liouville form itself does not change. 
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A general gauge transformation between two double-null gauges, U, V and U, V, 
reads 

U = X{U,o), V = Y{V,o), (71) 

where X and Y are suitable functions; the inverse transformation can be written as 

U = X{U,o), V = Y{V,o). (72) 

The form fU + gV + hid^ transforms under the change (0) of variables as follows 
fU + gV + hiY = fil +~gV + KY , 

where 

f = fx^u. g = ~gYv, 

and 

hi = hi + fX^i + gYi . 

For the function R{U, V, o) = -R(f/, V, o), we obtain 


R u — R t'tX, 




R.v = R i/Y, 


v^,v 


and 


The transformation of the logarithm is 

XuR' 


R.i — R.i + hijjXi + RyYi . 


In (-Q = In (- 

Collecting all terms, we obtain 


F = — 
dU 

G = — 
dV 


Hi = 


A 


R\ x,u 
'R\ f^u 
'X,u 


YyV 
R^X 


= In (-F) + In 


X 


,u 


Y 


V , 


uu 


+ FX, 


,u , 


R^ 


In 


Y 


,v 


4 Xy 

■,Ui Yyi 


4 

fX 

~4 


X 


,u 


Y 


(73) 

(74) 

+ H, + FX, + GY, . (75) 


y 


The first terms on the right-hand sides of Eqs. (^)-(^) represent a divergence; 
according to the result of the previous section, these terms can be thrown away, and 
we obtain finally: 

IT U'V R^Hyu 

F = FXy - -—^ , (76) 


G = GY 


' 4 Xy 
Yyv 


y 


4 Y 


y 


Hi = Hi + FXi + GY,-^{^ 

" 4 \Xy 


Y 


yi 


Y 


(77) 

(78) 


,v 
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These equations yield the transformation of the functions F, G and Hi that solve 
Eqs. (|^)-(|57[) for the gauge f/, V to the functions F, G, and Hi that solve them for 
the gauge U, V. 


3.3 Integration of equations ( |55|) — ([57| ) 

The transformation equations (|76D-([78D can be used to calculate the solutions F, G, 
and Hi from some well-known solutions F, G and Hi. In fact, one can always choose 
the gauge U, V in such a way that the right-hand sides of Eqs. (^ ) - (|T^) simplify 
and become trivially solvable. Let us, for instance, choose the gauge as follows. 

Let the volume term to be transformed correspond to the flat spacetime. Then we 
can choose the coordinates U and V as the retarded and advanced time coordinates 
for Minkowski spacetime and 


i = l , 


R 


-U + V 
2 


Eqs. now read 


Fy — Gjj = 0 , 

^i,U ~ ~ 0 ’ 

Hi Y — G^i = 0 , 


and we guess a solution to be 


F = G = H, = {) 


(79) 


for all i. 

Let the volume term correspond to the (a/3)-quadrant (see I, Sec. 2.1) of the 
Schwarzschild spacetime with the mass parameter M(o) (the mass parameter is 
Dirac’s observable and so it is, in general, a function of the chosen complete system of 
those observables). Then we can choose t/" and to be the double-null Eddington- 
Finkelstein coordinates dehned by Eqs. (1) and (2) of I so that Eqs. (3), (4) and (5) 
of I hold: 


A 

R 


2 M(o) 

R 


a(3 


2M{o)k 


a/3 exp ■ 


aU^ + (3Vf^\' 

4M(o) )\ ’ 


(80) 

(81) 


where n is dehned by Eq. (5) of Ref. |^. On substituting Eqs. (|80|) and (|8T|) for A 
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and R into Eqs. (|^-(|^, we obtain 


Again, we guess easily 


Fy Gjj — 0 , 

Ku-F, = 

Av-A. = 


F = 0, 

G = 0 , 


Hi 


aU“ + pVf ,, 

- - —M 


(82) 

(83) 

(84) 

(85) 


The double-null Eddington-Finkelstein gauge may be simple, but it is singular 
(t/" = aoo at the future and = —jSoo at the past horizon). The solution 
(P^) and (|5B[ ) diverges at both horizons. This singularity can, however, be removed 
by subtracting a suitable hh-term. To show this, let us transform to a regular gauge, 
for example to the Kruskal coordinates U and V: 


Then 


U = — exp 


/ aU^\ 
V 4M ) 


V = exp 



X(F, o) = -4aM(o) In(-F) , 
Y(F, o) = 4j3M(o) ln(l/) . 


An easy calculation using Eqs. ([7^) - ([78|) yields 


F 




—2MM j In 



These functions are singular at the horizons (F = 0 or V = 0). However, if we 
choose 

U-(C/,V',o) = -Af2(o)|n(^^^ _ 

then the equivalent solution dehned by 


Reg ■— F -h W^u , Greg G + Wy , 


Hireg Hi + W^i 


19 







is regular everywhere. We obtain (using some properties of the function k, see 


Eq. (12) and (51) of §) 

2exp(R/2M) ’ 

(86) 


^ M{R + 2M) „ 

“ 2exp(R/2M)^ ’ 

(87) 


Hireg 0 , 

(88) 

where 

R = 2Mk{-UV) . 



The solutions ([79D and (|83|)-(|85|) or (|86|)-(^ together with the formulae (|76D- 
( 0 ) can help us to calculate the functions F, G and Hi for the two-shell system. 
For the single shell, the resulting formulae of |0] can easily be reproduced by the 
new method. 
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